Collective interference of composite two-fermion bosons 
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The composite character of two-fermion bosons manifests itself in the interference of many com- 
posites as a deviation from the ideal bosonic behavior. A state of many composite bosons can be 
represented as a superposition of different numbers of perfect bosons and fermions, which allows us 
to provide the full Hong-Ou-Mandel-like counting statistics of interfering composites. Our theory 
quantitatively relates the deviation from the ideal bosonic interference pattern to the entanglement 
of the fermions within a single composite boson. 
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The quantum statistics of bosons is most apparent 
in correlation functions and counting statistics. Char- 
acteristic bosonic signatures are encountered for ther- 
mal states, which feature the Hanbury Brown and Twiss 
effect PHI], as well as in meticulously prepared Fock- 
states [3H6], which exhibit Hong-Ou-Mandel-like (HOM) 
interference. Deviations from the ideal bosonic pattern 
in HOM setups are often caused by inaccuracies in the 
preparation of Fock-states and in the alignment of the 
setup, which induce partial distinguishability between 
the particles @ 1 p. Another source for devia- 
tions from perfect bosonic behavior has received only 
little attention, limited to mixed states [12j HB]: Since 
most bosons are composites ("cobosons") made of an 
even number of fermions, reminiscences of underlying 
fermionic behavior are expected in many-coboson inter- 
ference. In analogy to partially distinguishable particles 
[21 HI], one can intuitively anticipate that the many- 
coboson wave-function partially behaves in a fermionic 
way, with impact on the resulting counting statistics. 

Here, we investigate such compositeness effects in 
HOM interferometry of cobosons. The ideal bosonic in- 
terference pattern is jeopardized by the Pauli principle 
that acts on the underlying fermions, an effect that be- 
comes relevant when the constituents populate only a 
small set of single-fermion states. The effective number of 
single-fermion states can be related to the entanglement 
between the fermions, via the Schmidt decomposition. 
Not only does entanglement thus guarantee the irrele- 
vance of the Pauli-principle for coboson states, but it also 
constitutes the very many-body coherence property that 
ensures that many-coboson interference matches the ideal 
bosonic pattern [3J [3] . The many-coboson wavefunction 
can be described as a superposition of different numbers 
of perfect bosons and fermions, with weights that are 
determined by the Schmidt coefficients. Using that in- 
tuitive representation, we compute the exact counting 
statistics in many-coboson interference and thus provide 
direct experimental observables for compositeness. Prop- 



erties of the collective wave- function of the fermionic con- 
stituents can thus be extracted from coboson interference 
signals, while in the limit of truly many particles, partic- 
ularly simple forms for the interference pattern emerge. 

The bottomline of our discussion, the observable com- 
petition of fermions for single-particle states, is a rather 
general phenomenon that is not restricted to any partic- 
ular physical system. To render our analysis of many- 
coboson interference tangible, however, we focus on an 
intcrfcrometric setup that can be realized with trapped 
ultracold atoms 1141. 
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FIG. 1: Setup for the interference of engineered cobosons. 
N\ (N2) strongly bound bi- fermions are prepared in the up- 
per (lower) lattice at J v <C Jh, such that each bi-fermion is 
governed by the local energies Ej and the tunneling rate Jh- 
The barrier between the lattices is then ramped down, such 
that J v ~S> Jh and vertical tunnelling takes place. The total 
number of bi-fermions in the upper and lower lattice is then 
counted. 

We consider strongly bound bi-fermion pairs that are 
trapped in a two-dimensional potential landscape with 
different horizontal and vertical coupling rates |15j . as 
depicted in Fig. [I] which is described by the Hamiltonian 
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where d q ■ = a q jb q ■ creates a bi-fermion consisting of an 
a- and a 6-type fermion in the jth site of the upper or 
lower lattice (q — 1, 2); Jh {J v ) is the effective tunneling 
strength along (between) the lattices, and tj defines a 
local energy landscape |14j . We assume that, initially, 
Jh S> J v , and mult i- cob oson states are prepared in the 
horizontally extended lattice q by [THl E] 
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A coboson is thus a horizontally dclocalized bi-fermion, 
and the S coefficients Xj are then the Schmidt coefficients 
of the two-fermion state. 

The distribution A is conveniently characterized by its 
moments 



M(m) = ^A™, 
j'=i 



(3) 



where normalization implies M(l) = f and M(2) = P 
is the purity of either reduced single-fermion state. We 
consider an initial state of N\ cobosons in the upper and 
N 2 cobosons in the lower lattice 1141 . 



which witness the probability to find exactly n\ (77,2) bi- 
fermions in the first (second) lattice. This procedure, 
however, is computationally expensive and does not offer 
an intuitive physical picture. By exploiting the symme- 
try properties of the state Q , one can show [2] that the 
behavior of cobosons is imitated exactly by a superposi- 
tion of states with a different number of perfect bosons 
and fermions, in analogy to partially distinguishable par- 
ticles [SJ [TT] . When the distribution of the bi-fermions 
along the lattices is neglected, exhibits the same total 
counting statistics in the two lattices as the state 
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where g q (f q •) creates a boson (j-type fermion) in the 
lattice q. The weight of the component with p pairs 
of fermionically behaving bi-fermions depends on the 
Schmidt coefficients and reads [14] 
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where we assume iVi > N%, and xn is the coboson nor- 
malization factor [T6ll20| . a symmetric polynomial [21] 
given by \n = H({1, . . . , 1}), with 
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To assess the behavior of the cobosons, we let the bi- 
fermions tunnel vertically between the two lattices by set- 
ting J v >• Jh and letting the system evolve for a time of 
the order 1/J V . Thus, beam-splitter-like dynamics cou- 
ples the two lattices, while tunneling processes within 
the lattices, induced by Jh, can be neglected on this time- 
scale. The Schmidt modes j are therefore left unchanged. 
Time-evolution until t implements a beam-splitter with 
reflectivity R = cos 2 (tJ v /2). In principle, the counting 
statistics of bi-fermions in the two lattices can be ob- 
tained by integrating the dynamics induced by Q for 
the initial state |\I / ) given in Eq. Q and taking the ex- 
pectation values of the counting operators 
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Combinatorially speaking, w p is the probability that, 
given two groups of Ni and N 2 objects with properties 
distributed according to A, and assuming that all ob- 
jects in either group carry different properties, one finds 
p pairs of objects with the same property when the two 
groups are merged. In the present context, w p denotes 
the population of the state components in which the 
Pauli principle affects p pairs of bi-fermions. The term 
|0(O)) thus describes perfect bosonic behavior, its weight 

wo = Xn 1 +n 2 /(^n 1 Xn 2 ) can be bound via the purity P 
and the particle numbers iVi , N 2 [13 [2D] : 

(L — Ni)\(L — N 2 )l 
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where L = [p] . 

We can now derive the counting statistics of cobosons 
after time-evolution until t — -k/2/ J v , which corresponds 
to a balanced beam-splitter with R = T = 1/2. The 
probability Pt t( m ) to find m cobosons in the upper lat- 
tice is the sum of the resulting probabilities from the 
different contributions in Q, 



Ptot(m) 



N 2 
p=0 



P(m,p), 



(11) 



where P(m, p) is the probability to find m particles of any 
species in the upper lattice, given the state \4>{p)) defined 
and the beam-splitter reflectivity R — 1/2 [T4"] . 
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FIG. 2: (color online) Counting statistics for the coboson- 
state |$) with JVi = N2 = 2, and of its components with 
different numbers of bosons and fermions \cj>{p)), p = 0, f , 2. 
Dark blue circles represent bosonically behaving bi-fermions, 
light orange symbols stand for fermionically behaving bi- 
fermions. The total counting statistics P to t (m) is the weighted 
sum ( f f I over the different components of the wave-function. 
While |(/>(0)) exhibits perfect bosonic behavior, \<j>(jp > 1)) are 
partially fermionic, which leaves a signature in the counting 
statistics. Here, R — 1/2 and Ai = • • • = A4 = 1/4, such that 
wo = W2 = 1/6, wi = 2/3. 



The simplest case is given by two interfering cobosons 
(iVi = N 2 = 1), for which we find wq = P and w\ = 1— P: 



P to t(l) = P, Ptot (0) = Ptot (2) = 



1 - P 



(12) 



For P — > 1, the Pauli principle dominates and one al- 
ways finds one particle in each lattice. In contrast to the 
interference of unbound boson pairs that can break up 
dynamically [22] , a perfect bosonic dip emerges here in 
the limit of vanishing purity, P — > 0. 

Higher-order power sums M(m) with m > 3 become 
relevant when more than two cobosons interfere. For 
example, the interference of N% = 1 with Ni cobosons 
reflects the normalization ratio Xn+i/xn [ISHI81 12"31 124) : 

Pot(m) = felp( m ,0) + ( 1 - M±±\ P( TO ,1).(13) 

x Nl V x Nl J 

In general, the balance between all the weights 
wq,...,wn 2 governs the counting statistics. Since the 
weights w p depend on power-sums M(m) up to or- 
der Ni + N 2 , the characteristics of the distribution A 
can be established through interference signals. For 
Ni = N2 = 2, we illustrate the decomposition ^ in 
Fig. [2] The ideal boson interference pattern P(m, 0) 
is jeopardized by the finite purity P = 1/4, the con- 
tributions of the single fermion-pair and double fermion- 
pair part in the wave-function lead to the altered signal 

Ptot(™-)- 

Distributions with the same purity P may have differ- 
ent higher-order power sums M(m), with consequently 



distinct counting statistics. Keeping P constant, the 
counting statistics is extremized by two particular distri- 



butions: the upper bound in ( 10 1 is saturated by peaked 
distribution A (p) with A^ p) >X 2 p) = ■■■ = \f\ in the 
limit S — > 00; the lower bound is saturated by the uni- 
form distribution A<") with < = ■■■ = X { ^ [1/p] , 
for fractional purities P = 1/L [20] . The counting statis- 
tics for iVi = N 2 = 6 is shown in Fig. [3| The weights 
vj^ pS> of the uniform (peaked) distributions differ consid- 
erably (see lower panel), which is reflected by the count- 
ing statistics (upper panel, note that P(m) = P(12 — m) 
due to symmetry). Only one Schmidt coefficient in the 

peaked distribution is finite in the limit S — >• 00, thus 

(p) (p) 

only the weights w and w\ are non-vanishing: the 
interference patterns of 12 and of 10 bosons take turns. 
Instead, all weights Wq" - <6 alternate for the uniform dis- 
tribution. Kinks emerge at fractional values of P, when 
a new non-vanishing Schmidt coefficient emerges. For 
P — > 1/6, fully fermionic behavior is attained, and one 
always finds six cobosons in each lattice. 
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FIG. 3: (color online) Upper panel: Counting statistics 
Ptot(m-) as a function of the purity, for the uniform (it) 
(right-hand part) and peaked (p) (left-hand part) distribu- 
tions \V*' P \ Lower panel: Corresponding weights «/ p /") of 
the coboson wavefunction given in 0. We set Ni = N 2 = 6, 
R — 1/2. The counting statistics is perfectly bosonic for van- 
ishing purity, P — > 0, while cobosons behave as fermions for 
the uniform distribution and P = 1/6. The number of non- 
vanishing Schmidt-coefficients in the uniform distribution is 
L = [1/P], hence the weights lUj with I < N — L — 1 van- 
ish: There are at least Ni — L — 1 pairs of fermions, which 
results in the kinks in the weights. The binomial distribution 
corresponds to the statistics of distinguishable particles. 

The dependence of Ptat(k) on the power-sums M(m) 
can be used to infer the latter from measured counting 
statistics for different N\,N%. The purity P follows im- 
mediately for iVj. = N% = 1 via Eq. (12 1; in general 
M(m) is inferred by the counting statistics of a total of 
N\ + N 2 — m cobosons. Since higher-order power-sums 



4 



are constrained by Jensen's and Holder inequalities [2"5] . 
M(m- 1)^=5 < M(m) < M(m- 1)^, (14) 

bounds for higher-order M(m) become tighter with in- 
creasing knowledge of M(m), as depicted in Fig.|4j 

M(m) 1.000: 
pm/2 0.500 



A(") with S Schmidt coefficients (P =1/5), 



S 
P 



o 



T3 

U 

N 

"3 

s 

o 



0.100 
0.050 - 

0.010 " 
0.005 " 

0.001 '_ 




■ Original distribution v 
Constraints by M(5) 
Constraints by M(4) 
Constraints by M(3) 

■ Constraints by M(2)=P 



4 5 6 

Power-sum order m 



FIG. 4: (color online) Normalized power-sums and con- 
straints. The normalization to p m / 2 is chosen such that the 
upper bound is constant. A randomly chosen distribution 
A leads to a certain hierarchy of power-sums (black stars). 
The measurement of interference signals with Ni and N2 co- 
bosons reveals the power-sums up to order JVi + N2, which 
leads to the indicated constraints on higher-order M(m) with 
m > Ni + N2 + 1 (blue, orange, green and red symbols), 
according to Eq. (1141). 



When the exact counting statistics cannot be retrieved 
and many (N > 1000) cobosons are brought to interfer- 
ence, such as in the interference of BECs the gran- 
ular structure of the interference pattern becomes sec- 
ondary. The impact of imperfect bosonic behavior can 
then be incorporated into a macroscopic wavefunction 
approach [3], i.e. the number of particles is treated as 
the amplitude of a single-particle wavefunction. Fock- 
states are modeled by a random phase between the dif- 
ferent components of the wavefunction. When the frac- 
tions Ij — Nj/(Ni + N 2 ) of ideal bosons are prepared in 
the two lattices, the particle fraction / in the upper lat- 
tice after beam-splitter dynamics obeys the probability 
distribution [3] 



Pmwf(P h,h) — 



1 



■Ky/mThh -(I- Rh 



for ARTIJ2 > {I - Rh - TI 2 ) 2 , while it vanishes oth- 
erwise. For cobosons, a finite fraction of fermions needs 
to be accounted for in each lattice. The probability dis- 
tribution for the particle fraction / then becomes 



V{I) 



dl f W(I f ) Vmwf (I - If, h-I f ,h- If) ■ 



where W(lf) is the probability distribution for the frac- 
tion of fermions If in each lattice. For the uniform state 
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The continuous limit W'"'(I/) is obtained for 
A4 + iV2 =: TV — ^ 00, when N\, N 2 ,p and S are scaled 
linearly with N: 

W^(I f ) = Im^ (N ■ w^ =If . N) ) =5 (I f - phh) ,(16) 

and the total number of bi-fermions per Schmidt mode is 
constant, p = N/S. Since the number of bi-fermions in 
either lattice is limited by S, it holds < p < l//i < 2. 
The fraction of perfect fermions is thus exactly the frac- 
tion of expected pairs of bi-fermions in the same Schmidt- 
mode, pl\l 2 , which gives 

V (U) (I) = Vmwf (I - phh\h{l - ph), h(l - ph)) ■ 

The width W of this distribution is closely related to the 
fraction of fermions, 



W = ly/BThhd - Ph)(l ~ Ph), 



(17) 



and becomes narrower with increasing number of bi- 
fermions per Schmidt mode, p. In principle, this 
may jeopardize Fock-state-interferometry with non- 
elementary particles such as neutral atoms, since the 
width of the intensity distribution is used to infer a small 
phase (which translates here to a reflectivity R). 

Trapped ultracold atoms typically feature very small 
electron-state purities of the order of 10~ 13 such 
that atom interferometers are not sensitive to the com- 
positeness of the atoms. With attractively interacting 
fermionic atoms in tunable external potentials , 
the transition between fully bosonic (P — > 0) and fully 
fermionic (P — > 1) behavior may be implemented ex- 
perimentally by varying the size of the available singlc- 
fcrmion space and observing the resulting interference 
pattern when bi-fermions are brought to interference |14) . 

In conclusion, even though two fermions may be ar- 
bitrarily strongly bound to a coboson with no appar- 
ent substructure, deviations from ideal bosonic behav- 
ior can be observable in many-coboson interference. Not 
the binding energy, but the entanglement between the 
fermions is observable on the level of the cobosons. 
The superposition Q allows to understand the partially 
fermionic behavior of cobosons, and ultimately leads to 



simple expressions for the interference of BEC (17). The 



methods that we have exposed can be extended immedi- 
ately to larger numbers of sublattices and to more com- 
plex interference scenarios j3]- 

Cobosons always constitute indistinguishable particles; 
two cobosons in the two lattices share the same distri- 
bution of Schmidt coefficients A. The impact of partial 
distinguishability and the effects of compositeness can 



actually be discriminated in the experiment: While par- 
tially distinguishable particles can be described as a su- 
perposition of perfect bosons and distinguishable parti- 
cles [21 [TT] , cobosons exhibit the behavior of a superpo- 
sition of bosons and fcrmions, which naturally leads to 
differing interference patterns in the two cases (see also 
the binomial distribution in Fig. [3j which is attained for 
distinguishable particles). 

The role of entanglement for bosonic behavior is 
twofold: It circumvents the Pauli principle for composite 
bosons [TMISl l20l 124] . and it maintains many-particle 
coherence. Quantum correlations between the fermions 
are necessary for the bosonic exchange symmetry in the 
relevant parts of the wave-function that allows the rep- 
resentation in Eq. Q . If mixed states of bi- fermions are 
prepared instead of entangled states, the exchange sym- 
metry and the encountered bosonic behavior break down 
- even though the combinatorial argument that relates 
to the number of accessible states remains valid. The 
visibility of correlation signals of, e.g. large molecules, is 
thus not only affected by the mixedness of the molecules 
at finite temperatures, but also by the consequent loss of 
many-particle coherence. 
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SUPPLEMENTAL MATERIAL 
Physical model, Hamiltonian Eq. (1), and preparation of the state of Eq. (4) 



For a tangible model of composite bosons, we consider fermions of two distinguishable species, a and 6, which 
interact attractively via a contact interaction U, and which are prepared in two weakly coupled one-dimensional 
lattices, as depicted in Fig. 1 of the main text. 

The single-particle tunnelling rates between the wells in horizontal (vertical) direction are denoted by Jo_h (Jo,v), 
such that the Hamiltonian on the level of the individual fermions reads 



2 S ^ 1 S 

h = - u J2J2 4 AA.Aj + 2 EE e > (4,Aj + KjKj 
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Jo,v 
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3 a %3 



+ h \ . h2 . +h .C.) , 



3 = 1 



where q = 1, 2 denotes the two one-dimensional sublattices, and j = 1, . . . , S the potential wells along each lattice. 
The attractive interaction U is always strong, U 3> Jq V i Jo h- Therefore, two fermions of the two different species are 
always bound to a bi-fcrmionic particle described by 



d) ■ = a* -P ■ 
q,3 q,3 q,3> 



which fulfills the algebraic properties of a hardcore-boson operator pQ , 

p ^ q or j : [<^ i3 -,4,» 



0. 



(ii) 



(iii) 



under the assumption that the numbers of fermions of each species coincide. 

Due to the strong attractive interaction between the fermions, they can only tunnel as a pair through off-resonant 
processes, with rates 



Jh = 



Jf},h 

~1T 



Jy - - 



T 2 

J 0,v 



(iv) 



for the interlattice and intralattice directions, respectively. We thus recover the effective Hamiltonian given by Eq. (1) 
in the main text: 



H = 



9=1 3 = 1 



3 = 1 



(v) 



9=1 3 = 1 



While the main purpose of our analysis is an investigation of compositencss for the collective interference exhibited 
by the state Eq. (4) in the main text, it is worthwhile to indicate a procedure by which such a state with many 
cobosons in the same state may be produced as the result of an experimental protocol. 

We may start with an extended lattice with a total of SNi sites (remember N\ > N2 by assumption), such that 
site kS (with k = 1, . . . , Ni — 1) is coupled to site kS + 1 with a strength J s -C Jh, and eus+j = e j f° r k= 1, . . . , N\. 
We exploit the exact mapping of hardcore-bosons to fermions in one dimension to obtain the ground-state of JVi and 
N2 bi-fermions in the first and second extended lattice as a direct product of the lowest Ni and N2 single-particle 
states, respectively, 



|GS(JVi,JV a )) 



Ni /SNi 

HE 

3=1 \2=1 



^3\l d \,l 



N 2 /SNi N 

n e 

j=i \i=i j 



|0), 



(vi) 



where the matrix Uji contains the coefficients 1 = 1... SNi of the jth single-particle eigenfunction in one sublattice. 
After preparing the ground state with N\ and N2 particles in each sublattice, we project away the component of the 
many-body wavefunction in which particles are present in the sites S + 1 . . . SN\ of either lattice. Consequently, a 
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particle initially prepared in the jth eigenstate in lattice q is projected - with a finite probability - onto the state 
created by 

s 

E^ d L- ( vii ) 

{=1 

Our above choice of the potential landscape with J s <C Jh ensures that the JVi energetically lowest single-particle 
wave-functions have a node around / = S and that the coefficients ujjj with 1 < I < S are very similar for the first 
Ni eigenfunctions, 

2 < j < Ni,l <l < S : u) jt i « cjij. (viii) 

After projecting out components that populate the auxiliary sites j = S + 1 . . . SN\, the many-particle-state is 
consequently very close to a state of N q -io\d population of the co-boson state 

s s _ 

4 = a E w i.i4,j = E v^i«?,i' ( ix ) 
i=i i=i 

where a ensures normalization and we thus set A; = lacji^l 2 . Note that the Pauli principle ensures that unlike for 
hardcore bosons the preparation of fermions into this state is impossible: One then always finds at least one particle 
in the sites S + 1 . . . SNi, and the projection always fails. 

In the main text, we take the initial state (4) as given, and thus consider coboson operators for the two lattices 
q = 1,2 as given by Eq. ( |ixj ). By appropriately modelling the local energies tj, the resulting distribution of Schmidt 
coefficients A; can be modelled to a very wide extent. 



Behavior of cobosons under beam-splitter dynamics, Eqs. (7,8) 



We would like to describe the counting statistics of the state of Eq. (4) in the main text in an efficient manner. 
For that purpose, we insert the definition of the coboson creation operator, Eq. (nxh , into the initial state, Eq. (4): 



JV 2 



l<kj,lj<S r ffi 



E 



n 



l.fer, 



iV 2 



n x^'A 



2./„ 



(x) 



All indices k m (l n ) appertain to the upper (lower) lattice. It may occur that k m = l n for some m, n, i.e. two bi-fcrmions 
can occupy the same well j in the two different lattices. The sum Q can be written in terms with a given number of 
pairs of indices p that fulfill k rn = l n . There can be between none and N2 of such pairs (remember N\ > N2): In the 
former case, k m ^ for all m, n; in the latter, k m = l m for all m < N% (disregarding permutation of indices). The 
state initial state thus becomes 



N 2 

|*)=£|$(p)) 

p=0 



(xi) 



where 



I*(P)> ■= 



N 2 
P 



'(N 2 ,P) 

E 



n 



i,fe„ 



N 2 

n 

n=l 



P J XJVi • N^. ■ xn 2 ■ N 2 l 
and the sum Y^^ 2 '^ runs over all indices k m , l n (1 < m < Ni, 1 < n < N 2 ) that fulfill 



|0) 



(xii) 



i ^ j 
1 < m < p 
p < m < N\,p < n < N 2 



ki 7^ fcj 1 7^ lj 5 
k — I 

kjn ~/~ In • 



(xiii) 
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Reordering indices and setting N := Ni + N 2 — 2p, we can rewrite the sum (xii) as 



I*(P)> = 



N 2 \ ( Ni 




V J XN t ■ Nil ■ xn 2 ■ N 2 l 



l<lj<S Vl,m:njtl m 

E E E 









N 




-m- P 










n \J x °(o) 


















_ rn—1 





N 2 -p 

ih 



'2,a(Ni-p+n) 



form ionic 



bosonic 



(xiv) 



where the indices r\ 
By inspecting 



replace the k m>p and l m>p and S 



? } denotes the permutations of the 



we can now infer the time-evolution of each |$(p)) component of the many-coboson state ( pdj ), 

contains p pairs 



as it is induced by the Hamiltonian (1) in the here-considered limit J/j <C J„: Each summand in 
of bi-fermions that occupy the same well in the two lattices - they take into account the summands with k 
(1 < m < p) of Eq. 



Due to the Pauli principle, these bi-fermions cannot tunnel and thus behave fermionic. The 

(which correspond to the indices n,m > p in Eq. 



other bi-fermions described by the indices r±, ...,r^ (which correspond to the indices n,m > p in Eq. (xii)) always 
occupy different wells, such that the Pauli principle does not apply and tunneling is possible. One such summand of 
is depicted in Fig. |lj Since the state of the bi-fermions that can tunnel is fully symmetric under the exchange 



of any two bi-fermions between the lattices (note the sum over all permutations of the lattice indices in (xiv)), it is 
manifestly bosonic, as also illustrated in Fig. In] 




ri 



>2 



r?, 



FIG. I: One summand of Eq. ( |xiv[ ), with iVi = 4, A?2 = 3: Two pairs of bi-fermions (p = 2) are located in the same wells h and 
h, and cannot tunnel. All other bi-fermions can tunnel, and interfere with other terms in the sum (see Fig. [lip . 



©=0 iO=G 
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Q=® 
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(a) 




(b) 



FIG. II: Emergence of Hong-Ou-Mandel-like bunching for bi-fermions. Each component of the wave-function in ( |xiv[ | , corre- 
sponding to a permutation 01, interferes with another component 02 in which the bi-fermions swap the wells they occupy, 
(a) Destructive interference: The two processes (both bi-fermions remain in the same lattice, or both bi-fermions tunnel) lead 
to the same final state with one particle in each lattice. The right-hand-side process, however, acquires a phase of i 2 = — 1 
due to tunneling, such that the two processes interfere destructively, (b) Constructive interference: The final state with both 
bi-fermions in the same lattice is fed by two processes with one tunnelling event, i.e. both paths acquire the same phase and 
constructive interference takes place. 



The total state \&(p)) thus exhibits the same counting statistics as a state \<t>{p)) of p pairs of distinct fermions and 
N% + N2 — 2p perfect bosons (Eq. (9) in the main text): 



n 



Nq—P 



L = W(^-P)! 



n 



J l,3 



10} , 



(xv) 



9 



where g| (fqj) creates a boson (j-type fermion) in the sublattice q while the actual location along the lattice is 
omitted (remember that the number of bi-fermions in each lattice is counted, independently of the location of the 
bi-fermions along the lattice). Since no interference between different p occurs, the initial state |\E') (Eq. (4)) behaves 
like a superposition of states \4>{p)) with different numbers of ideal fermions, in analogy to partially distinguishable 
particles [2]: 

N 2 

=^2 ^/w^\(f>(p)) , (xvi) 

p=0 

where 

Wp = (*(p)|*(p)) = ( Nl \( N f) -^—n({2,...,2, 1 1 !). (xvii) 

\P J \ P J XN 1 XN 2 v v ' " v ' 

p Nx+N 2 -2p 

is the weight of the component with p pairs of fermionically behaving bi-fermions, and Q({xi, . . .xn}) is given by 
Eq. (5) in the main text. 



Evaluation of counting statistics, Eq. (11) 



The counting statistics that is exhibited by the substitute state Eq. (xvi I, can be inferred by inserting the 
single-particle time-evolution for the creation operators for each state \(j>(p)), 

g\ -> iVRgl + Vfgt_ q 

where R = cos (J v t/2) and T = sin 2 {J v t/2) are the reflection and transmission coefficients of the beam-splitter 
dynamics. 

The probability P(m,p) to find m particles in the upper lattice can then be inferred by taking the overlap of the 
state \4>{p)) after time-evolution with the Fock-state of m — p bosons and p fermions in the first mode, 



m—p f + ^ N—m—p 

n /iV3,j io> ■ 



WW 



y/(m - p)\(N - m - p)\ j=1 
The evaluation of the overlap can be done following the methods presented in Refs. [2HJ]. 
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